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In this work, we have considered modified Chaplygin gas with variable G and A. The trivial 
solution describes decelerating phase to accelerating phase of the universe. The non-static with 
constant equation of state describes the inflationary solution. For static universe, G and A must be 
forms arbitrary and for static universe with constant equation of state G and A should be constant. 

PACS numbers: 

I. INTRODUCTION 

The Einstein field equation has two parameters - the gravitational constant G and the cosmological constant 
A. The Newtonian constant of gravitation G plays the role of a coupling constant between geometry and 
matter in the Einstein field equations. In an evolving Universe, it appears natural to look at this "constant" 
as a function of time. Numerous suggestions based on different arguments have been proposed in the past 
few decades in which G varies with time [1]. Dirac [2] proposed a theory with variable G motivated by the 
occurrence of large numbers discovered by Weyl, Eddington and Dirac himself. Many other extensions of 
Einstein's theory with time dependent G have also been proposed in order to achieve a possible unification of 
gravitation and elementary particle physics or to incorporate Mach's principle in general relativity [3]. 

From the point of view of incorporating particle physics into Einstein's theory of gravitation, the simplest 
approach is to interpret the cosmological constant A in terms of quantum mechanics and the Physics of vacuum 
[4]. The A term has also been interpreted in terms of the Higgs scalar field [5]. Linde [6] proposed that A 
is a function of temperature and related it to the process of broken symmetries. Gaspirini [7] in this regard 
argues that A can also be interpreted as a measure of temperature of a vacuum which should decrease like 
the radiation temperature with cosmic acceleration. By considering the conservation of the energy-momentum 
tensor of matter and vacuum take together, many authors have invoked the idea of a decreasing vacuum 
energy and hence a varying cosmological constant A with cosmic expansion in the frame work of Einstein's theory. 

A as a function of time has also been considered by several authors in various variable G theories in different 
contexts [8]. Investigating the distance dependence of gravity under very general conditions, Wilkins [9] found 
that the gravity field at a distance r from a point mass has two components: one varying as the other as 
r (Hookian field). The latter component is identifiable with the weak field limit of the A term in Einstein's 
equation. His analysis allows one to consider both the gravity fields - the Hookian field, coupled to A and the 
Newtonian one coupled to G - on an equal footing. With this in view, several authors [10, 11] have proposed 
linking the variation of G with that of A in the framework of general relativity. This approach preserves 
conservation of the energy-momentum tensor of matter and leaves the form of the Einstein field equations 
unchanged. Though this approach is non-covariant, it is worth studying because it may be a limit of some 
higher dimensional fully covariant theory [10, 12]. 

Recent observations of the luminosity of type la Supernovae indicate [13, 14] an accelerated expansion of 
the Universe and lead to the search for a new type of matter which violates the strong energy condition i.e., 
p -I- 3p < 0. The matter consent responsible for such a condition to be satisfied at a certain stage of evolution 
of the universe is referred to as a dark energy. There are different candidates to play the role of the dark 
energy. The type of dark energy represented by a scalar field is often called Quintessence. The simplest 
candidate for dark energy is Cosmological Constant A. In particular one can try another type of dark energy, 
the so-called Chaplygin gas which obeys an equation of state like [15] p = —B/p, {B > 0), where p and p are 
respectively the pressure and energy density. Subsequently the above equation was generalized to the form [16] 
p = —B/p'"-,0 < n <1. There are some works on modified Chapligin Gas obeying equation of state [17, 18] 
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p = Ap-B/p-,{A>Q) (1) 
At all stages it shows a mixture. This is described from radiation era to ACDM model. 

II. EINSTEIN'S FIELD EQUATIONS 

We consider the homogeneous and isotropic space-time given by FRW metric 



1 — fcr^ 



87rG(t)p A(t) k 



and 



Prom (5) and (6) we have 



a^'^ ^' G 8ttG 
We assume the law of conservation of energy {T^- = 0) giving 



(2) 



ds^ =dt^ -a^{t) 

where k (= 0, ±1) is the curvature parameter. 
The energy-momentum tensor for perfect fluid is 

Tij = {p + p)uiUj + pgij (3) 
where p and p are energy density and isotropic pressure respectively and c = 1. 
The Einstein field equations with variable G and A is given by 

Rij - \ R9ij - mg,, = -87rG(t)r,, (4) 

i.e., we have two equations as 

a2 



+ ^-7^ (5) 



a2 3 3 a2 



- = -^(p + 3p) + Ai (6) 



X G A „ 
p + 3-{p + p)+p - + — = (7) 



p + 3^{p+p) = (8) 

From (7) and (8) we have 

A = -87rGp (9) 

This implies G > or < according as A < or > i.e., G increases or decreases according to whether A 

decreases or increases. 



Now for modified Chaplygin gas (1) we obtain 
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Now from equation of state (1), field equation (5) and conservation equation (8) we have 



(10) 



^=9( 1 + A- 



B 



3 3p a?p 



Differentiating w.r.t t, we obtain 

f2p 3/32 \ 3B(n + l)/33 

^ V ^ ) p"+3 (l + ^ - ^) 



= 3p 1 + A - 



B 

r,n+l 
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2k 



(11) 



(12) 



Now we consider p ^ (tlie case p = i.e., /? = constant is discussed later). Since the above differential 
equation of p is highly non-linear, so it cannot be solve analytically. Now for simplicity of calculation, we choose 
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2p 



3B(n + l)p2 
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= 



From equation (13) we find the trivial solution of A i.e.. 



(13) 



(14) 



A = 



(1 + 2>A)k _ 2kB 3(i+A)(i+n)-2 



(1 + A)a^ C{\ + Af 
Also from equation (14) we obtain (after manipulation) 



(15) 



3(1 + A) 
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(16) 



For fc > 0, A decreases with t upto certain stage of the evolution of the universe and for fc < 0, A in- 
creases with t after certain stage of the evolution of the universe (see figure 1) and for /c = 0, we must have A = 0. 



From (9) and (15) we have 



(l+3A)fc /cB{3(l+A)(l+«)-2} 3{(i+A)(l+n)-l} 



B C 



(17) 



Note that here Go depends on the value of k i.e., k = Q implies Go = 0. Prom figure 2 we see that G increases 
with the evolution of the universe for fc = 1. 

For early universe i.e., for a w 0, we get 



47rG w -5^ 1 + 47rGo , (Go > 0) 

(1 + A)G— 



(18) 



4 



20 



- 




k=l 









/ ; 






/ ' 






/ 






/ 






/ : 






/ 






/ ' 

/ 



FiK.l 




Fig.2 




1 2 3 4 5 

t 

Fig.3 

Fig. 1, 2 and 3 show the variations of A, G and g against t iox A = 1/3, B — \, C = 1, n = 1/2. 



i.e., G — ^ Go as a — 0. 

Also for late universe i.e., for a « oo, we have 



4^G « 7,-^^ ( a3(i+^)(i+")-2 + 4nGo (19) 



C(l + A)2 V B 
i.e., G — cx) as a — )■ oo. 

The deceleration parameter has the expression 

ad 4ttG{p + 3p) - a 



SttGp + a - 3fca-2 



(20) 



From figure 3 we see that q decreases from +1 to —1 for A — 1/3, k = ±1. This implies the universe has 
early deceleration and late acceleration. For fc = 0, we have A = and G — 0. So the field equations (5) and 
(6) yield a — constant, which is static solution. 



The above discussions are valid for p ^ i.e., p ^ constant. Now we will discuss in the case p = i.e., p = 
constant = po (say). 
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Now equation (8) reduces to 



ip + p) - = 



(21) 



So we have possibilities: p + p = Oord = 0. 



(i) d 7^ and p + p = 0: This implies p = —p = —po = Po (say). From equation (1) we have po = 
So from field equations (5) and (6), we get 



B \ 1+" 

1+^ . 



cosh(-\/Ci t) 



Prom equation (9) we have 



A = Ao - SttG 



B 



1 + A 



(Ao = 3Ci) 



(22) 



There arc no other equations, so A and G can not be calculated. We also see that in this case, G increases 
or decreases according as A decreases or increases. A and G are arbitrary functions of time in this case. 

The field equations (5) and (6) and equation (22), we have 



t = 



log 



HV^ a + VAoa^ - 3fc 



/A^ ao + \/Aoao - Sk 



I.e., 



a = ao cosh(y — t) 



3k . Ao 



(23) 



For fc = we have a = ao ev 3 * which is the inflationary solution. 



(ii) d = and p + p ^0: This implies a = constant = oo i.e., we have static universe. In this case, from fleld 
equations we have the values of G and A as 



G 



k 



Anal 



(1 + A)po 



Po 



— constant 



(24) 



and 



A = 



(1 + 3A)p^+^ - W 



= constant 



1 

where po ^ • 

(iii) d = and p + p = Q: This implies a = constant = ao i.e., we also have static universe and p = 

1 

'Po = Po (say). Prom equation (1) we have po = (1^) "^^ • From field equations, G and A satisfies: 

3k 

A + SnGpo = 

do 



(25) 



-P = 



(26) 
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and 



A + SnGpo = 



(27) 



which are consistent only for k — 0. This imphes 



A 



SttGpo = SttG 



( 



l + A 



B 



) 



+1 



(28) 



The above relation (28) shows that G and A are arbitrary functions of time t. 



III. CONCLUDING REMARKS 



From the field equations and the conservation equation, an equation is obtained in p, a and A which suggests 
two trivial solutions. The trivial case with A is given in equation (15) leads to a model which starts from big 
bang with non-zero gravitational constant Go {k ^ 0) and has a positive constant deceleration parameter q = 1 
and leads to accelerating universe (g = —1) with infinite gravitational constant G (see figure 3). For fc = 
we have A = and = and a = constant, which is static solution. The model describes the evolution of 

the universe with early deceleration and late acceleration. For constant density (po 7^ (tTa) ^ ) universe, 
the cosmological constant A and gravitational constant G are arbitrary functions of time in which for static 
universe describes fc = 0. But for static universe with non-constant density, G and A are constant where 
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